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k> : Abstract 

^ Superheavy (M > 10^'^ GeV) particles produced during inflation may be the dark matter, inde- 

pendent of their interaction strength. Strongly interacting superheavy particles will be captured 
by the sun, and their annihilation in the center of the sun will produce a flux of energetic neutrinos 
that should be detectable by neutrino telescopes. Depending on the particle mass, event rates in 
a cubic-kilometer detector range from several per hour to several per year. The signature of the 
process is a predominance of tau neutrinos, with a relatively flat energy spectrum of events ranging 
from 50 GeV to many TeV, and with the mean energy of detected tau neutrinos about 3 TeV. 
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I. INTRODUCTION 

It is usually assumed that the interactions of dark matter and ordinary matter are weak, 
at least as small as ordinary weak interactions. It is also often assumed that the dark matter 
is a thermal relic of the big bang. If it is a thermal relic particle, then its mass should be less 
than 340 TeV, the unitarity limit ^. These considerations lead to the popular picture for 
relic dark matter of an an electrically neutral particle, without strong interactions, and with 
mass less than a couple of hundred TeV. In this paper we explore a path less traveled, and 
assume that the dark matter is a nonthermal relic, it interacts strongly with normal matter, 
and it is very massive. We show that the clean signature of this possibility is a detectable 
flux of energetic neutrinos from the sun 

It has long been appreciated that if the dark matter is massive, say larger than a few 
TeV, it will behave effectively as dissipationless dark matter regardless of whether it has 
strong or electromagnetic interactions 0, |^. However, since the upper limit to the mass of 
a thermal relic is a few hundred TeV, the window for very massive dark matter particles 
was thought to be not very wide. 

The recent development of scenarios for nonthermal production of dark matter has opened 
the window to the possibility that the dark matter might be supermassive, independent of 
its interaction strength 0, §, §. Of the many possibilities for producing supermassive 
dark matter, perhaps gravitational production is the most general @, §[. In this scenario, 
dark matter is produced by vacuum quantum fluctuations toward the end of inflation. The 
resulting particle density is independent of the interaction strength of the particle, which 
leads to the possibility that the dark matter may be electrically charged, strongly interacting, 
weakly interacting, or may have only gravitational interactions with normal matter. 

A particularly promising mass range for gravitational production of dark matter is the 
mass scale of the inflaton, about 10^^ GeV in chaotic inflation models. If the inflaton mass 
heralds a new mass scale, then it would be reasonable to imagine that there are other 
particles of similar mass. Furthermore, gravitational production of particles with a mass 
comparable to the inflaton mass naturally leaves behind a cosmologically interesting density 



of dark matter today 0]. The particle content may include exotic quarks or other strongly 
interacting particles in the spectrum of new particles . 

In this paper we will consider the case that the dark matter is strongly interacting and 
supermassive, a simpzilla. Although our calculations will not be sensitive to whether the 
SIMPZILLA is electrically charged, there are arguments that suggest that the simpzilla must 
be neutral P, Q. The possibility that the dark matter may be very massive and strongly 



interacting was recently discussed by Faraggi, Olive, and Pospelov |jTO|- In that paper they 
have a nice discussion of the particle physics motivations for the existence of a massive, 
stable, strongly interacting particle, and they point out that the sun and Earth may be the 
source of high-energy neutrinos from the annihilation of the particles. 

In the next section we will calculate the trapping rate and annihilation rate of simpzillas 
in the sun (and Earth) as a function of the simpzilla mass and interaction cross section. In 
Section III we will discuss the emergent spectrum of neutrinos from simpzilla annihilation 
in the center of the sun. In Section IV we will calculate the event rate in cubic-kilometer 
underwater or underice neutrino detectors. Finally, the last section contains our conclusions. 

A preview of our main conclusion is that we expect a detectable high-energy solar neutrino 
flux for much of the parameter range of interest. 

II. CAPTURE AND ANNIHILATION RATE 

Before launching into the details of the capture-rate calculation it is useful to make some 
extremely crude estimates. The first estimate is for the number of SIMPZILLAS that hit the 
sun. 

Assuming the simpzillas comprise the local dark matter density of 0.3 GeV cm~^, the 
local number density of simpzillas is 3 x 10"^^ M^^^cm"^, where M12 is the simpzilla 
mass in units of lO^^GeV. Assuming a typical velocity of 240 km s~^ for the simpzillas, 
the local flux is approximately 6 x lO^^M^^^cm^^s^^. The surface area of the sun is about 
6 X lO^^cm^, and the product of the surface area and the flux, 4 x 10^^Mi2^s^^, is a crude 
estimate of the rate of simpzillas hitting the sun. 



Now consider trapping of simpzillas in the sun. Assuming the simpzilla impacts the 
sun with the solar escape velocity of 600 km s~^, the kinetic energy of the simpzilla is 
about 10^Mi2GeV. Suppose the simpzilla scatters with nucleons with a cross section 
of 0" = 10^^'^cr_24 cm^ |ll|], and in every collision suffers an energy loss of mproton'y^/2 = 



2 X lO^^GeV. Through the center of the sun is about 2 x lO^^cm"^ of material, so there will 
be about 2 x 10^^a_24 colhsions with a total energy loss of 4 x 10^o"_24GeV. Since the initial 
kinetic energy is about 10^Mi2GeV, for masses much less than 10^^o"Z24GeV, simpzillas 
with average velocity hitting the sun will be trapped, but if the simpzilla mass is much 
larger than that, only the low- velocity tail of the phase-space distribution will be captured. 

These considerations determine the gross behavior of the dependence of the simpzilla 
trapping rate on the simpzilla mass and its interaction cross section. For a mass of about 
lO^^GeV, about 10^^ simpzillas hit the sun per second. For an interaction cross section of 
about 10~^'*cm^, most of those simpzillas are captured. 

Now we turn to the details of the capture calculation. The capture of dark matter particles 
in the sun and Earth has been studied in detail [1^, 0; here we adapt the considerations 



to the capture of simpzillas by the sun. The capture rate of dark matter by the sun 
is given in terms of an integral over f{u), the phase-space density of dark matter in the 
solar neighborhood normalized such that J^ f{u)du gives the number density of particles 
at R, some sufficiently large radius where the gravitational pull of the sun is negligible. The 
capture rate is given by 

1 



c 



A-kW 



f{u) udu dsw? 6 



p{e,u), (1) 



where 6 is the angle between the velocity of the particle and the normal to the surface at 
R, and hence the expression inside the square bracket is the contribution to the inward flux 
at radius R from particles at velocity u and angle 6. Finally, P{6, u) gives the probability 
that particles with the given velocity and direction at some large radius R will be captured 
by the sun. 

Writing the angular momentum per unit mass as J = Ru sin 6, and defining Wq = u'^ + Vq 



where Vq is the escape velocity at the surface of the sun, we can rewrite the above as 

poo pRqWq rj 72 

Tc = Ti f{u)udu -^PUu), (2) 

Jo Jo u^ 

where Rq is the solar radius, and P{J,u) now gives the capture probability for a particle 
with angular momentum J and a velocity at infinity of u. 

The regime in which dark matter capture is often considered is that of weakly interacting 
particles, where the optical depth of the sun is much less than unity, i.e., njycri^Q <^ 1 where 
n^- is the number density of nuclei in the sun, and a is the dark-matter — nucleon cross 
section. In this case P{J,u) is given by |1T2| : 

\Jo u^ J Jo Jo u^ ^i _ ^jjrwY \ 4^2; \^^ wy 

where w now stands for u"^ + v"^ where v is the escape velocity at radius r, and fi = mx/rriisi 
where tjix is the mass of the dark matter particle (in our case, the simpzilla) and ttitv is 
the average nucleon mass. The first term in square brackets is a normalizing factor, the term 
2nN<jdr/ Jl — (J/rwY gives the probability of collision, and the terms in the second line give 
the probability that the particle gets scattered into a bound orbit. The step function Q{x) 
equals 1 or depending on whether its argument is positive or not. It is simple to show that 
a collision produces a fractional energy change AE/E that is uniformly distributed between 
and 4//i where we have taken the large fi limit (otherwise, 4//i is replaced by Afi/{fi + 1)^). 
The above assumes that the dark matter particle typically suffers at most one collision in 
its passage through the sun. We are here interested in the opposite regime where tinctRq 3> 
1. For each collision, the fractional energy change is of the order of l//x. Therefore, after 
unctRq collisions in the sun, the simpzilla would get captured if Vq ^ (fg + u'^){l — 
1/yu)"^'^^®, implying all particles with M < m^. get captured where m^, = f0[l/(l — l//i)"'^'^'^° — 



1]-*^/^. For capture by the sun, it adequate to approximate m* by u^, = Vq/ J fi/{nNaRQ) — 1 
if fi/rixaRQ > 1 or m* = 00 if h/unctRq < 1. In other words, we approximate P in Eq. (H) 
as a step function. 

Using the Maxwell-Boltzmann phase space distribution function for the simpzillas 

f{u) = A{nx/J'Tcul^){u/uth)'^exp{—u'^/u'^Yi) [0] with px = 0.3(GeV/mx)cm~'^, we obtain 
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FIG. 1: The capture rate by the sun (solid curves) and Earth (dashed curves) for different 
values of the siMPZlLLA-nucleon cross section. The three curves for the sun and for Earth 
correspond to, from top to bottom, a = 10^^^, 10"^"^, and 10~^^ cm^. 

from Eq. (H) the capture rate for SIMPZILLAS. The capture rate has two forms, depending 
on the efficiency of energy loss in the sun. This is parametrized by q, defined as 



(1 = 



/^ 



20 



(4) 



unctRq "VlQi^GeVyV a /VTxlOiOcmy V ^0 
If g < 1, the SIMPZILLA will be efficient in losing energy in its passage through the sun, and 
the capture rate is 



Tc = 10'' s-'il + y 



lO^^GeV 



Uth 



mx J V240kms-V V7 x IQi^cm 



R, 



'0 



(5) 



where 

^~ ■ V600kms-V V240kms~V ' ^^ 

On the other hand, if g > 1, only low velocity simpzillas are captured and the capture rate 
is 



(7) 
where 



X 



(8) 



The capture rate as a function of simpzilla mass and cross section is shown in Fig. ^ 

The captured simpzillas settle into the core of the sun on a time scale determined by 
their drift velocity, tdrift ~ "'"/^drift where f drift can be estimated by balancing gravity with the 
viscous drag due to scattering with nucleons [|I^]: Gpr^mx/r'^ ~ crraArfdrift"^Af(^/'"^Af)"'^'^^- 
Here, p and T is the typical mass density and temperature of the sun. The timescale is very 
short, of order 100 s {a/lQ-^^ cm'^){10^^ GeV/mx)if/10'' K)^/^. 

The subsequent evolution of the collection of simpzillas at the core can be divided into 
two stages. The first stage is when A^, the total number of simpzillas, is less than Nsg, the 
critical number necessary for the SIMPZILLAS to become self gravitating. In this first stage 
the SIMPZILLAS are supported by the thermal pressure of the surrounding plasma and the 
SIMPZILLA profile follows an isothermal distribution given by pxi'f') oc exp(— r^/2r^) where 



r, = (STQ/AnGmxpQY^^ = 5000 cm {W^ GeV /mxY^'^{TQ/lO'^ Ky^'^{150 gcm-^ p^. 



1 1/2 



and Tq and p© are the core temperature and density respectively. The critical number 
is defined by 47rr^p0/3 ~ Nsg^^x, and is equal to 



■l0^2GeV\ 5/2 /^^x 3/2 / i50gen i-3\ ^'^ 
mx ) \Wk) I p0 



AT^^^lO-p-rr^ _|_ p^l^r:^ . (9) 



To determine if N ever exceeds NsG) "we have to determine Neq-, the number of simpzillas 
in the sun when there is an equilibrium between annihilation and capture. If the size is 
given by r* above, then Neq is found by equating the annihilation rate and the capture 
rate: 2(o"yif )(A''£;Q)^/(47rr^/3) ~ Vc- Using the annihilation cross-section (cr^f) < l/{rn\v), 



we obtain 

The above means that unless mx is significantly smaller than lO^^GeV, no equilibrium is 
reached for A^ < Nsg- With the capture rate given above, N can reach Nsg on a time-scale 
much shorter than the lifetime of the sun. The next stage is then set for the collapse of the 
SIMPZILLA collection, whose final state will be determined by a number of factors. 

First, the critical A/'chandra, beyond which the collection of simpzillas cannot be sup- 
ported by degeneracy pressure, is given by balancing Gm^-A^'chandra/''^ ~ '^^^^^cdandra/'''- 

7VcH„.,.~a'/^10-(^3^)"'. (11) 

where a is the ratio of electron number density to simpzilla number density. This is 
initially a large number, i.e., the mass density of simpzillas and the mass density of 
nuclei are initially comparable, which implies a ~ mx/rriN. It is unclear how much a 
would be reduced in the collapse process. It depends on how effective simpzillas are in 
dragging along protons. If a > 10^°(mx/10^^ GeV), A^grav is smaller than A^chandra, and so 
the collapse would result in a simpzilla collection supported by non-relativistic electron 
degeneracy pressure. If a ^ 10^°(mx/10^^ GeV), the configuration will collapse to a black 
hole, unless sufficient annihilation occurs along the way. 

Let us consider the question of whether annihilation would halt an otherwise catastrophic 
collapse. Following 0, the fractional change in A^ can be estimated by 

^^ - r,fdt3^!^. (12) 



A^ J Anr^t) 

Whether sufficient annihilation occurs or not depends on whether the integral is dominated 
by small r (late times) or large r (early times). The longer the configuration spends at small 
radii, or in other words, the slower the acceleration, the better the chance for annihilation 
to work against collapse. The relation dr/dt ~ (ro/r)^(ro/tdrift) was used in Ref. 0, where 
tq is the initial radius and tdrift is the drift-time given before (the timescale for collapse is set 
by viscous drag). A perhaps more reasonable limit to how fast a given shell can accelerate 



is given by free fall: dr/dt oc 1/r^/^. Integrating, one obtains: 



N 27r rl'\]''' ' 



(13) 



rf - 10-^ cm (^)'" ( irl X' C-^f' C-^^^) (14) 

^ \WV I 10^32 cni3 s-1 / VlOOs/ V rn / ^ ^ 



where we have taken the limit of rj -C tq. Clearly, if Vf were sufficiently small, AN/N ~ 1 
can always be achieved. The only thing one has to make sure is that the required Vf is larger 
than the Schwarzschild radius. Using ro ~ r* and A^ ~ Nsg, it can be verified that 

2/3 

/ Mrift \ ' / JU 

ro 

would do the job. 

After shedding a fair fraction of the SIMPZILLAS, the equilibrium configuration should in 
principle be one where support is provided by non- relativist ic electron degeneracy pressure: 
Gm\N/r ~ a^/^A^^/^/(2mer^) where rrie is the electron mass. The above, together with the 
requirement of the balance of annihilation and capture, gives 

Comparing the above with A"chandra shows that the final configuration is just barely stable, 
depending somewhat on the exact value of a and nix- If not stable, then the configuration 
goes through another cycle of collapse and eventual halt by annihilation. It is curious that 
this cycle might go on indefinitely, in which case each collapse would be accompanied by 
enhanced annihilation and therefore a mild neutrino outburst (using for example AA^ ~ 10^° 
and tdrift ~ 100 s gives an annihilation rate of 10^^ s~^, not overwhelmingly larger than the 
capture rate of ~ lO^^s"^). We will assume for the rest of this paper that the annihilation 
rate is given by capture rate, or more precisely, 2Va ~ ^c- 

Finally, there is the possibility that a drops to a sufficiently small value later on that 
electron degeneracy pressure is irrelevant, and the simpzillas are supported instead by 
their own degeneracy pressure (or its analog if it were a boson; see e.g., Ref. [T^). In this 
case the expression equivalent to Eq. (|TT]) is 

iVchandra ~ ^^^(lO^^ GeV/m^)^, (16) 



and the equilibrium number will be 

Neq - 10^(rc/10i' 8-^)^/^(10^2 GeV/mx)'/'(^;/10-^c)^/^ 

so that the final equilibrium configuration is clearly stable. 

For completeness, we give the capture rate of simpzillas by Earth: 



(17) 
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X 10^2^-1 



l+y - exp(-a;2)(l + y^ + x^ 



mx 



where now y = 0.04, and x is given by 



x = y 



1 



(l-/i 



-i\K. 



1/2 



N,on = 2 X 10" 



a 



fi= 10 



12 



mx 



(18) 



• (19) 



, 10-24 cmV "^ VlOi^GeV. 

The parameter A'^cou is the average number of collisions the simpzilla suffers. It may be 
identified with fi/q in the case of the sun |]l^. The expression for Tc is well approximated 
by the corresponding expressions in Eqs. (^ through (0) for y > 1, which applies for the 
case of the sun. For Earth we have to resort to this more complicated expression; however, 
it has simple limiting forms. For A^cou//^ ^ 1? one can use [cf. Eq. (^] 

lO^^GeVX 



c 



X W^s-\l+y' 



nix 



For Ncowl H ^ 1, the capture rate is well approximated by [cf. Eq. (|^)] 

.. . . . , ,,^12 



X 10^2 s-^ x^y^ ( — 



mx 



(20) 



(21) 



with X = y/\/q — 1, where g = 5 x 102(mx/10-'^2QgV)(10'24cm2/cr). For most, but not all, 
of the parameters of interest, it is this last regime that is relevant for capture by Earth, 
which gives F^ ~ 4 x 10''s"^(10^2QgV/mx)2(cr/10-24cm2). The total number of simpzillas 
captured in the lifetime of Earth, Ie ~ lO^^s, is then 4 x 10'^^{10^'^GeV/mx)'^{a/10-^'^cm^). 
We can compare this with the number of simpzillas in equilibrium using Eq. ([lOD , with 
the replacements Tq —>■ T^ c^ 5000 K, p© ^ p^ ~ 10 g cm~^, and t> ^ 2 x lO^^^c, obtaining 
Neq ~ 5 X 10^^(mx/GeV)^/2(p^/400 s^'^Y^'^. Therefore, the SIMPZILLAS captured in Earth 
would be in equilibrium unless the mass or cross section is significantly different from what 
is assumed above. In general, the annihilation rate is given by: 



2 X 10^"^ 
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mx 
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4 X 10% 
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lOi^s 



.(22) 



The capture rate for Earth is shown in Fig. [^ along with the capture rate for the sun. 



III. SIMPZILLA ANNIHILATION IN THE SUN 

Take a simple picture where the SIMPZILLA annihilates and produces two quarks or two 
gluons, each of energy mx where mx is the mass of the SIMPZILLA. The quarks and gluons 
then fragment into high multiplicity jets of hadrons and secondary decay products. 

Defining x = E/Ejct = E/rrix, one can take the fragmentation function for the total 

number of hadrons to be [|1^] 

^ = ax--'/\l-xf. (23) 

ax 

Here a is some constant that can be set by total energy considerations: 

1 = / x^T^dx = a-—, (24) 

JO ax 15 

so a = 15/16. 

The total number of hadrons produced in the fragmentation of the jet is 

"1 dN„ , 15 1 



N, 



H 






(25) 



When calculating the total multiplicity, the cutoff for the integral should be e = Aqcd/"^x- 
Using Aqcd = 0.1 GeV and defining M12 = mx/lO^^ GeV, 

Nn='4(^T = Qx^0'Mli\ (26) 

The final decay chain of the hadrons will contain all species of neutrinos. 

We will be interested in the number of heavy quarks, bottom and top. If all quarks 
were light, then all flavors would be produced equally. However, because of the mass of 
the heavy quark, the number must be found by using e = Mq/mx in Eq. (|26D rather than 
e = AqcD/mx- 

where Mq is the mass of a typical meson containing the heavy quark: 2 GeV for the 
charm, 5 GeV for bottom, and 175 GeV for top. Therefore, per annihilation, one expects 
Nch^rm/ Nh = 0.23, Nbottom/ Nh = 0.14, and Ntop/Nn = 0.024. This means that each 
annihilation into two jets produces 7.4 x lO^M^g light hadrons, 2.8 x lO^M^g charmed 
hadrons, 1.6 x 10^Mi2 bottom hadrons, and 2.8 x lO^M^l top hadrons. 
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FIG. 2: The spectrum of top hadrons produced by the fragmentation of quark and gluon 
jets from SIMPZILLA annihilation (dotted hne), and the v^ spectrum produced by decay of 
top (dashed hne). For the fragmentation spectrum, -Emin = 175 GeV and A^total = 2.8 x 
10^M]^2 psr annihilation. For the decay spectrum, we will only be concerned with neutrinos 
above -Emin = 50 GeV. The number of tau neutrinos above -Emm is ^Vtotal = V^'^M^^ P^r 
annihilation, and half that for the other neutrinos. 



It is also possible to estimate the spectrum of the heavy-quark hadrons: 



-En 



dN \ { E \-^l'^ 



{E>E, 



(28) 



A'^TOTAL dE 2 V-Em: 

where -Emin is the mass of the top quark, bottom quark, or charm quark. 

The resulting -E~^/^ fragmentation spectrum for top hadrons is shown by the dotted line 
labeled "fragmentation" in Fig. |^. The minimum energy is approximately the mass of the 
top quark, 175 GeV. 

We will assume that simpzilla annihilation occurs in a medium of density found in 
the center of the sun, p ~ 200gcm~'^, or n ~ lO^^cm"^. Using an interaction cross sec- 



tion of 10"^^ cm^, the hadronic interaction length is about 10~^cni. Light and charmed 
hadrons scatter many times before decay and the resultant neutrinos will have very low 
energy. The B lifetime is about 10"^^ s, so the decay length is 3 x 10^"^ {E/ Mb) cm, or 
6 X 10~'^(-E'/GeV) cm using Mb = 5 GeV. The ratio of the decay length to interaction length 
is Ld/Li = 0.6(-E'/GeV). So for E > 5 GeV the B will also scatter and lose energy before 
decay. It won't completely stop after one scattering. Most of the cross section is diffractive 
production of low-energy crud and there should still be a leading B, so B decay is a potential 
source of high-energy neutrinos. 

However, top hadrons are a promising source. The top lifetime is short, and almost 100% 
of the time decays as t — > Wb. The W lifetime is 3 x 10^^^(£'/80GeV) s, which results 
in a decay length of about 10^^^(£'/GeV)cm, which will be much less than the interaction 
length for even the most energetic tops. The W then decays with a branching ratio of 1/3 
as W —^ lui (equally into tu^-, jiu^ and eVe)- Therefore it is reasonable to assume the top 
quark will produce energetic neutrinos before losing energy. 

The spectrum of neutrinos produced in the chain t ^ W ^ v is straightforward to 
calculate. A convenient analytic fit to the spectrum is given by 

'^^M , ^ + ^^»^ (29) 



where A/" is a normalization factor. This spectrum is shown in Fig. ^. As expected, at 
energies larger than the top mass the E"^/"^ spectrum is recovered. 

We will use the spectrum in Fig. |] as the spectrum of neutrinos produced by simpzilla 
annihilation. Since there are about 2.8 x 10^M]^2 tops produced per annihilation, and 10% 
of them make TVr followed by r decay including a z/,-, the total yield of z/t-'s per annihilation is 
5.6 X lQ'^Mi2 ■ Top (as well as r) decay also produces /i + z/^ and e + z/g, but the electrons are 
absorbed and the muons are stopped before decay, so the yield of high-energy z/^ and u^, per 
annihilation is about 3.8 x lO^M^g ■ Only about 20% of the neutrinos will be produced with 
energy above 50 GeV, so the emission rate of tau neutrinos in the core is about 10^Mi2 ^a 
with a spectrum above 50 GeV shown in the figure. This, of course, is the emission rate in 
the core. We now turn to the propagation of the neutrinos through the sun. 



IV. THE EMERGENT SPECTRUM OF NEUTRINOS 

The total neutrino emission rate above 50 GeV from the core of the sun is /core = 
10^Mi2 Ta for Vri and half that for z/^ and z/g. Here T a = rc/2 is the simpzilla annihilation 
rate. The core emission rate spectrum above 50 GeV for each neutrino species is [see Eq. 

(H)] 

(^] = 333 GeV^/^ ^ 



dE E, 



mm 



""^^^ Q{E-E^,^), (30) 

^[E + Mt] [{E + MtY - Mf] [{E + Mwf - M^] 

where -Emin = 50 GeV and the theta-function vanishes for negative argument and is unity 

for positive argument. Of course df /dE has been normalized such that 

™'° \ / core 

But since the sun is opaque to energetic neutrinos, the emergent emission rate spectrum 
is not the same as the core emission rate spectrum. The emission rate spectrum of neutrinos 
that emerge unscathed is 

fdf\ 



V / unscattered 



i^) ^^P (-^(E) f n{r) dvY (32) 

V / core 



where n{r) is the radial dependence of the number density of the sun and acc{E) is the 
energy-dependent charged-current cross section given in Table |. Also given in Table | is the 
neutral-current cross section. 

The electrons and muons produced by charged-current interactions are rapidly thermal- 
ized, so the effect of charged-current interactions is effectively to remove u^ and Uf, neutrinos 
above a transparency energy E,^ where cr{Ef^) J n{r) dr becomes unity. Using 

n{r) = 1.4 X lO2^exp(-r/O.li?0) cm"^ (33) 

for the density profile and adopting the cross section for E < 10^ GeV, the transparency 
energy is E,^ = 150 GeV, and the emission spectrum of z/^ and u^, emergent from the sun is 

(S) ri^) '^''-'"'^- *''' 

\ / emergent \ / core 



TABLE I: The energy-dependent cross sections used to calculate the flux of neutrinos from 



the sun ^9 
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FIG. 3: The emergent emission rate spectrum of neutrinos from the sun. For i/^, /„ 
10^Mi2 r^, and for Un and Ug, half that value. In all cases, -Emin = 50 GeV. 



This emission rate spectrum is shown in Fig. ^ 

The situation is different for tau neutrinos. The hfetime of the tau produced in a charged- 
current scattering is so short that it decays before significant energy loss. Since the decay of 
the tau includes a Vr, the effect of the scattering and subsequent decay is simply to reduce 
the energy of the incident Vr to about 20% of its incident value. Therefore, an incident tau 
neutrino above the transparency energy will continually suffer tau production and decay 
interactions, but will not be removed from the flux of high-energy neutrinos. The process 
will continue until the tau-neutrino energy has been degraded to the transparency energy 
or below, then the tau neutrino will escape. 

This process has been considered by Halzen and Saltzberg ||2^ for very energetic neutri- 
nos propagating through Earth. They found that the emergent spectrum of neutrinos was 
very well described by a lognormal distribution centered on the transparency energy with a 
dispersion of 0.49 decades in energy. We will assume that the flux of scattered tau neutrinos 
above the transparency energy emerges in a lognormal distribution peaked at -E^ = 150 GeV 
with a dispersion of 0.49 decades in energy. 

The emergent emission spectrum of tau neutrinos has two contributions. The first con- 
tribution is the unscattered emission spectrum, given by Eq. (|32|). The second contribution 
is the fraction of the original emission above the transparency energy which will emerge as 
a lognormal distribution centered on the transparency energy; 

{%] = 1^ ^ exp [- (log E-\ogE^f /2a'] , (35) 

\"'^ / scattered -^m™ 

where a = 0.49 and JF is found by demanding that the integral of Eq. ( ^5|) results in the 
total number of neutrinos above 150 GeV that are scattered. (If the initial energy of the 
neutrino is below 150 GeV, a scattering will produce a neutrino below -Emin-) The result is 
J^ = 4.6 X 10-2_ 

So for tau neutrinos, the emergent neutrino spectrum is 

(4)_ =©_*©._■ 

This spectrum is also shown in Fig. & 



V. THE EVENT RATE 

In the last section we calculated the emission rate and emission rate spectrum of neutrinos 
from the sun. In this section we will calculate the event rate in a suppositious underice or 



underwater neutrino detector of approximate size of a cubic kilometer |21|. We will only 
consider the rate for "contained events," where the neutrino converts inside the volume of 
the detector. Including "uncontained events" will not significantly alter our results because 
the muon range at the relevant energy range here is comparable to a kilometer. We will 
assume that the efficiency of detection is a step function: zero below 50 GeV and unity 
above 50 GeV. 

The first step in calculating the event rate is the simple step of converting the emission 
rate spectrum calculated in the last section to a fiux spectrum arriving at Earth. Since the 
sun-Earth distance is D = 1.5 x 10^ km, the fiux spectrum is 

\ / emergent \ / emergent 

The mean-free-path of neutrinos is much larger than the size of the detector, so the 
fraction of the incident neutrinos of energy E that convert inside the detector is riicc o'cciE) L 
where nice is the number density of the target and L = 1 km is the size of the detector. The 
event-rate spectrum is given simply by 

^ = — [riiee acciE) L]AQ{E-50 GeV), (38) 

where A is the area of the detector, assumed to be 1 km^, and the G-function represents the 
detector efficiency. 

Since nice <^cciE) L = 4 x 10^^^{E/1 GeV), we find for the event-rate spectrum 

i - " >< 1«-" G^ (i) ,e(^-»GeV). (39) 

\ / emergent 

The total event rate is the integral of the event-rate spectrum. 

For muon or electron neutrinos, the total event rate can be found using Eq. ( p4D for the 
emergent emission spectrum, with result 

Rtotal{i^^, i^e) = 1.6 X IQ-^^mU'Ta. (40) 
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FIG. 4: The event rate in a cubic-kilometer underice/underwater detector for different 
values of the SIMPZILLA-nucleon cross section: from top to bottom the curves are for a = 
10~^^, 10~^^, and 10~^^ cm^. The upper solid curves are for tau neutrinos and the lower 
dashed curves are for muon and electron neutrinos. For comparison, the 90% CL. upper 
limits on muon fluxes of nonatmospheric origin in the direction of Earth's core or the sun is 
approximately 10^ km~^ year~^ ||22|| . 

This event rate is shown in Fig. |, and the event-rate spectrum is shown in Fig. |^. 

For tau neutrinos, the total event rate is found using Eq. (p6|) for the emergent emission 
spectrum. For tau neutrinos the total event rate is 



i?TOTAL(/^r) = 1.1 X W^Mti'TA. 



20a^1/2i 



(41) 



This event rate is also shown in Fig. ^, along with the event-rate spectrum in Fig. ^ 

The contribution of scattered tau neutrinos dominates the event rate. While the u^^ and 
z/e events would be peaked toward the lower energy of the detector and drop rapidly, the Ur 
events would have a relatively flat spectrum extending from the lower limit of the detector 
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FIG. 5: The. spectrum of events for different neutrino species. The total event rate 
-Rtotal = 1.1 X 10^^°Mi2 Ta is the sum of the electron-neutrino, muon-neutrino, and 
tau-neutrino rates. Also shown are the two contributions to the tau-neutrino events, from 
unscattered neutrinos and from scattered neutrinos [see Eq. (p^)]. 

out to about 1000 GeV. 

The mean energies of the detected neutrinos are 



{E)= 197 GeV {u^,Ue) 
(E) = 3454 GeV (z/,). 



(42) 



We close this section by remarking on the possibility of a detectable event rate from 
annihilation of simpzillas captured by Earth. The ratio of the event rate for neutrinos of 
solar origin to the event rate for neutrinos of terrestrial origin is shown in Fig. |^. For small 
mass or large interaction cross section the signal from the center of Earth my be larger than 
the solar signal. For a < 10~^^ cm^ the solar signal will dominate for simpzilla masses 
larger than about 10^ GeV, while for a < 10^^^ cm^, the solar signal dominates for the 
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FIG. 6: The ratio of the event rates from solar and terrestrial neutrinos originating from 
SIMPZILLA annihilation. The ratio is a function of SIMPZILLA mass mx and scattering cross 
section a through the dependence of the capture rates on mx and a. 

entire range of SIMPZILLA mass considered here. 
VI. CONCLUSIONS 



If the local dark matter is very massive and strongly interacting, it should collect in the 
sun and Earth in sufficient numbers that equilibrium will be maintained between the capture 
rate and the annihilation rate. 

Annihilation or decay of very massive particles into hadronic channels in the solar core 
or at the center of Earth leads to the production of high-energy neutrinos. While electron 
and muon neutrinos above 150 GeV are mostly absorbed, energetic tau neutrinos will be 
emitted, and will be the signature of energetic jet fragmentation at the center of the sun. 

For most of the range of parameter space, the event rate expected in kilometer-scale 
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FIG. 7: The shaded region above the jagged hne is excluded by a variety of considerations 
as discussed in 0. Shghtly stronger, mo del- dependent constraints can be found in Refs. 
P, ^. The shaded region above the straight line would result in more than 10 events per 
year in a cubic-kilometer underice or underwater detector and should easily be able to be 
excluded. 



neutrino detectors will be well within detection limits. Such detectors should be able to 
exclude (or confirm!) the possibility of simpzillas as dark matter. 

In this paper we have only considered production of neutrinos through top quark produc- 
tion and decay. Another potential source of tau neutrinos is bottom quark production and 
decay. While this may increase the emission rate of tau neutrinos, the spectrum is expected 
to be the same. 

It is important to note that there is essentially no background. For example consider the 
background from cosmic-ray produced neutrinos. To produce a 1 TeV neutrino, a center-of- 



momentum energy of around ^/s = J2mprotonE = 10 TeV is required, where E is the cosmic 
ray energy. Thus, a threshold energy of Eth = 5 x lO^^eV is required to produce a 1 TeV 



neutrino. The cosmic-ray flux at high energies is approximately 



2 



^' ^ = 10^'-' -^^ (43) 

ah s sr m^ 

which when integrated to give the total flux of particles with energy greater than E gives 

F{E > Eth) = TT^ 10^^-^ -^ = /-.^ , . (44) 

2Ej.fj s sr m^ yr deg km 

Since the sun subtends about a square degree, this sets the scale of the background. For 
much of parameter space, the signal should be well above the background. 

In Fig. ^ we present out results in the a vs. rrix plane and compare then with other limits. 
Clearly our considerations greatly extends the excluded region. 

Indirect detection of SIMPZILLAS through annihilation in the sun is complementary to 
the other idea for indirect detection: wimpzilla decay producing ultra-high energy cosmic 
rays H H. 

Finally, we comment on the possible role of neutrino oscillations. Neutrino oscillations 

will be important if 

Am2 lO^GeV L , , 

10-3eV^ E^ 2 X lO^cm ~ ' ^ ^ 

where L is the path length. Since 2 x lO^cm is Rq/?)^, if Am? is greater than or of the order of 
3 X 10~^eV^, oscillations will occur in the sun. If Am? is greater than or of order of 10~^eV^, 
then oscillations will be important during the neutrino's transit to Earth. Oscillations 
between v-r and v^ or z/g in the sun will decrease the z/,- emission rate, while oscillations in 
transit will change the flavor signature of the signal. 
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